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Validity of Transport Energy in Disordered Organic Semiconductors 
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A systematic study of the transport energy in disordered organic semiconductors based on variable 
range hopping theory has been presented here. The temperature, electric field, material disorder and 
carrier concentration dependent transport energy is extensively discussed. We demonstrate here, 
transport energy is not a general concept and invalid even in low electric field and concentration 
regime. 
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Understanding the charge transport mechanism in dis- 
ordered organic semiconductors such as conjugated and 
molecularly doped polymers, is of crucial importance to 
design and synthesize better materials. Physical trans- 
port phenomena in organic semiconductors are in gen- 
eral very complex and arise from energy and spatial dis- 
order, which invalidates the use of the traditional band 
transport theories based on the periodic distribution of 
atoms in inorganic crystals. Currently, charge transport 
in organic semiconductors is often described in terms of 
hopping theory A fundamental work in the model- 
ing of the hopping transport in organic semiconductors is 
the Gaussian disorder model proposed by Bassler 0, Q . 
Subsequent to this Monte Carlo simulation, a number 
of theoretical investigations concerning charge transport 
in amorphous organic systems utilizing the Gaussian ap- 
proach have been published and considerable progress in 
the analytical description of the problem has been made 
[i]-0] ■ An important simplification of the complex hop- 
ping transport mechanism was the introduction of the 
transport energy concept, This concept allows the com- 
plex hopping mechanism in the band tails to be inter- 
preted in terms of a multiple-trapping-and-release model, 
where the transport energy plays the role similar to the 
mobility edge in amorphous inorganic semiconductors 0- 
flij ]. Thermally stimulated luminescence carrier de- 
pendent mobility [lfl, |l5Hl7|, seebeck coefficient [18|, as 
well as injection phenomena[l9j| . have already been de- 
scribed utilizing this transport energy concept. Despite 
the concept of transport energy has widely been applied 
to describe the charge hopping transport of different or- 
ganic semiconductors, the validity of this concept is non- 
trivial. As transport energy is derived on the basis of zero 
electric field and low carrier concentration regime, and 
only the hopping upwards transport is concerned. How- 
ever, in the real disordered organic semiconductors de- 
vices, the electric field could reach 1 x 10 s V/m and carrier 



concentration could as high as 1 x 10 



20, 



, therefore, the 



applicability of the transport energy to real disordered 
organic semiconductors should be a matter of intensive 
research, where it is important to distinguish whether the 
organic semiconductor exist the mobility edge (transport 



energy) or whether this concept could be used to de- 
scribe the mobility in hopping system. In this letter, 
we show how the transport energy changes with temper- 
ature, electric field, carrier concentration and material 
disorder. Using this analysis, we conclude that the trans- 
port energy is not a general concept and invalid in real 
disordered organic semiconductors. 

Model. — Generally speaking, there are two ways to 
define he concept of transport energy in organic semi- 
conductors, both of them are based on the Miller- 
Abrahams expressions [20], the hopping transport takes 
place via tunneling between an initial states land a target 
statesj.The tunneling process is described as 



u exp (-it) = v a 



exp 2aRi 
(-2aRij). 



kT 



, Ei > Ej 
Ei < Ei 



(1) 

Here, Vq is the attempt-to-jump frequenct, Rij is the hop- 
ping distance, u is the hopping range 9],Ei and Ej are 
the eneries at site i and j, respectively, a is the inverse 
localized length and k is the Boltzmann constant. Then, 
if only the hopping upwards iAf is taken into account, the 
transport energy is defined as the finial energy E t that 
the hopping transport has maximum rate, so E t could be 
obtained by the equation as 12 | 

dEt 

This gives 
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g(E t 



E t 



9(E) dE 



1 



xkT 
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Where g (E) is the density of states (DOS). Here we can 
see that the hopping transport is limited by upward tran- 
sitions from filled states to empty states and the trans- 
port energy has been denned as the preferred energy to 
which the fastest upward transitions occur. 
In the other case, the transport energy is derived based 
on the average number of target sites n (Ei, u) as 



n (Ei,u) 



u/2a 



RijdRij 



Ei+kT{u-2a) 



g(E)dE. (4) 
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If we disregard the hopping downwards and choose the 
starting energy Ei as — oo, the transport energy could be 
obtained by setting n (Ei , u) = 1 and the result is 9] 



/ * g(E)(E t -E) 3 dE=-(akT) 3 



(5) 



This definition hints that the transport energy is the site 
energy to which the hopping upwards need the least en- 
ergy in energy space. 

When there exists an electric field F, the electric field 
will lower the Coulomb barrier, which leads to the reduc- 
tion of the thermal activation energies, and the hopping 
range with normalized energy can therefore be rewritten 
as I2lll22l 



2a (1 + /?) Rij + ej -e i ,e j >e i - cos f 
(2aRij ) . tj < €{ — (3 cos 9 



(6) 



Where (3 = Fe/2akT and and 9 is the angle between Rij 
and the electric field ranging from to pi. For a site with 
energy a in the hopping space, the most probable hop 
for a carrier on this site is to an empty site at a range u, 
where it needs the minimum energy. Conduction is the 
result of a long sequence of hops through this hopping 
space. Then, following the method in [21(, we derive the 
number of empty sites enclosed by the constant range u, 
as 



n(ei,u,F) = 



1 r r 

— =■ / d9sm9 dr2nr 2 
8« 3 Jo Jo 



i+e<-r(l+0< 



dc 



xg(e) [l-f( e ,e F )} .(7) 

Here / (e, ep) is the Fermi Dirac distribution and 1 — 
/ (e, 6f) is the probability that the finial site is empty, 
the Fermi energy ep is calculated by the condition 



dc 



9(e) 



1 + exp (e — ep) 



(8) 



Here c is the carrier concentration. By change the inte- 
gration variable, equation (7) will be in the form of 



n(e,u,F) 

€i~\-U 



2tt 



deg (e) 



3 x 8a 3 
(u + e. 



dr 



(1+j8t)- 



deg (e) 



(9) 



The first term on the right-hand side of equation (9) 
gives the number of shallower states, and the second one 
describes the number of target states which are deeper 
than the starting site . According to the variable range 
hopping theory |(| 3, 21 1, at a given field and temper- 
ature, almost every starting localized state has only one 
well distinguished nearest empty target hopping neigh- 
bor, i.e. another localized state that is characterized by 



the minimum value of the hopping range, this range could 
be obtained by solving the equation 



n (ei,F,u) = 1. 



(10) 



This equation established the basis of our model, the 
target site energy for every hopping process can be well 
evaluated from the euqation (10). 

Validity of Arkhipov transport energy. — To check this 
model validity, we disregard the downwards hopping as 
well, equation (9) then reads as 



1 = 



2vr 



3 x 8a 3 



dr I deg (e) { - + & ' ~ % (11) 
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Using the developed model equation (11), we now pro- 
ceed to calculate the relation between starting energy and 
finial energy in real disordered organic semiconductors. 
We take the Gaussian form g(e) = J^-^ exp 

density of states in the full manuscript, where N t is the 
number of states per unit volume and a = a^/kT indi- 
cates the width of the DOS. N t = 1 x 10 28 m~ 3 is taken 
in the full manuscript, a typical value for the relevant 
organic semiconductors. It is instructive to calculate this 
target energy as a function of the initial energy for dif- 
ferent carrier concentration, corresponding to different 
fermi-level, the results are displayed in the insert of Fig. 
1. A clear observation is that, in this situation, the trans- 
port energy does exist for deeper starting energy but 
increases with the carrier concentration. Field depen- 
dent transport energy is plotted in the Fig.l, the same as 
concentration dependent transport energy, field does not 
change transport energy for deeper energies but decreases 
for field higher than 1 x 10 7 V/m, this field strength is ac- 
tually low field for most organic devices. The reason for 
the decrease of transport energy is that, the electric field 
can change the energy difference between the finial and 
target sites, and thereby, assist carrier jumps along the 
field direction, so the target energy will decrease on the 
contrary. For the very low carrier concentration, the field 
and concentration dependent transport energy could be 
derived as 



n(e,u,F) 



2tt 



24a 3 (l-/3 2 r Jet-p 



g (e) (e, + u - e) de 



= — 2 / ' 9 (e) (et " e) 3 de (12) 

24a 3 (1 - /3 2 ) 2 J -00 ' V ' 

Certainly, this is not a general concept as well. 
We then consider the effect of the lattice disorder, which 
is caused by the random molecular packing in organic 
semiconductors j23|. The relation between transport en- 
ergy and materials disorder a is shown in Fig. 2 (a). 
An important result is that, the transport energy is in- 
valid for higher energies, but these energies make sense 
for the charge transport in real organic semiconductors. 
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It is well known, in the low carrier concentration regime, 
the hopping usually jumps from the so-called equilibrium 
energy = - (a/kT) 2 = -25fcT for a/kT = 5; In the 
high carrier concentration regime, the carrier jumps from 
Fermi energy {e F = —20kgT here). An obvious feature 
appears here is the transport energy does not exist for en- 
ergy higher than — 18/cT in the case of a/kT = 5. In or- 
ganic semiconductors, a — 5kT correspond to a = 0.12ev 
at room temperature, a typical value for organic semicon- 
ductors 24 |. Therefore, the transport energy has no gen- 
erality and is invalid for real organic semiconductor sys- 
tem. Please note that the field used here is 1 x 10 5 V/m, 
which is low field regime; the fermi energy chosen here 
is also corresponding to the low concentration regime. 
The other parameters values are also typical ones for or- 
ganic semiconductors such as a — lnm. To investigate 
the reason for this phenomena, we plot the ratio between 
downwards and upwards hopping in Fig. 2(b), one can 
see clearly that, the downwards hopping has little effect 
on the charge transport characteristics and should not 
account for this result. One possible explanation is, for 
the small disorder Gsussian DOS, since the carriers have 
few nearest neighbor sites to choose and it is reasonable 
to hop to the same energy; But for the wider Gaussian 
DOS, carriers has more near neighbor sites and more ac- 
tive, hence the target energy will be random accordin g to 
the calculation. Therefore, the approximation in [fl Il2j 
that the the transport energy in zero field is independent 
the starting energy is not correct in real organic semicon- 
ductors. 
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FIG. 2: (a) The disorder dependent transport energy for dif- 
ferent starting energy, (b) The ratio of hopping-downwards 
site number and hopping upwards number 
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FIG. 1: The computed transport energy as function of the 
starting energy for different electric field . The inset shows 
the carrier concentration (Fermi level) dependent transport 
energy. 



The question arises now on how much difference be- 
tween the calculated mobility using Akhipov method and 
the work here. The comparison is shown in Fig. 3. It is 
clearly seen that, in the low temperature regime, Akhipov 
model will underestimate the mobility but overestimate 
the mobility in high temperature. And this trend is even 



more pronounced with carrier concentration increasing. 

Validity of Baranovskii transport energy model. — 
According to the definition of Baranovskii transport en- 
ergy, the hopping upwards rate has maximum value for 
every jumping, hence, we derivative equation (6) as 



20(1 + 0)^ + 1 = 0, 
J ^ = 2a(l + /3)^=0, 
1 = AnRl f° sin OdO deg (e) [1 - / (e, e F )] . 

(13) 

The Rij that is the function of 6 and €j, could be ob- 
tained as the equation (11). Connecting Gaussian DOS 
and equation(12), we numerically calculate the transport 
energy for this definition. The concentration dependent 
transport energy is shown in Fig. 4. The deviation of tar- 
get energy from the transport energy is more dramatic, 
even at the starting energy e; = 24fcsT, the transport 
energy is obvious invalid. This holds for the field depen- 
dent transport energy (the insert figure), too. We also 
find, the same conclusion is obtained by introducing the 
percolation parameter [l2| in equation (13) . 

In conclusion, the validity of the transport energy in 
disordered organic semiconductors has been investigated 
intensively. The results shows, neither Baranovskii nor 
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FIG. 3: The Comparison between the calculated mobility in 
this work (F = 5 x W 4 V/cm) and Akhipov mobility model 
for different temperatures. 
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FIG. 4: The computed Baranovskii transport energy as func- 
tion of the starting energy for different electric field . The 
inset shows the electric field dependent Baranovskii transport 
energy. 



Arkhipov definition for the transport energy is valid 
in real organic semiconductors system, even in the low 
field regime, the transport energy lost the universality. 
This issue was not adequately addressed for in earlier 
evaluations of the transport energy based on the hopping 
rates. Concomitantly, the use of the previously obtained 
expressions for transport energy in calculations of the 
carrier transport parameters would lead to incorrect 
results for the concentration dependencies of these 
parameters [3, 25 - 27 1 . It should be mentioned that 



these calculations are describing only the first release 
step. These oscillating jumps of the released carriers 
may jump back to their initial state do not contribute 
effectively to the charge transport [9j. Finally, we 
mention that the description presented here should also 
be applicable to describe the temperature, electric field 
and carrier concentration dependent charge trans por t 
properties, for example mobility, as we have done in 
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